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Abstract 



The dynamics of magnetic bubble solitons in a two-dimensional isotropic antiferromagnetic spin lattice 
is studied, in the case where the exchange integral J{x,y) is position dependent. In the near continuum 
regime, this system is described by the relativistic 0(3) sigma model on a spacetime with a spatially 
inhomogeneous metric, determined by J. The geodesic approximation is used to describe low energy 
soliton dynamics in this system: n-soliton motion is approximated by geodesic motion in the moduli space 
M„ of static n-solitons, equipped with the metric 7. Explicit formulae for 7 for various natural choices 
of J{x, y) are obtained. From these it is shown that single soliton trajectories experience refraction, with 
analogous to the refractive index, and that this refraction effect allows the construction of simple 
bubble lenses and bubble guides. The case where J has a disk inhomogeneity (taking the value J+ outside 
a disk, and J_ < J+ inside) is considered in detail. It is argued that, for sufficiently large J+/J„ this 
type of antiferromagnet supports approximate quasibreathers: two or more coincident bubbles confined 
within the disk which spin internally while their shape undergoes periodic oscillations with a generically 
incommensurate period. 

1 Introduction 

The purpose of this paper is to analyze the dynamics of bubble-like topological solitons arising in a two- 
dimensional classical Heisenberg antiferromagnet where the exchange interaction is isotropic, but not homo- 
geneous. At each site («, j) € in a square lattice, one has a classical spin vector G M'^, with |Sy p ~ s^, 
evolving according to the law 



where r is time and Jij are positive constants. This is the usual planar isotropic Heisenberg antiferromagnet, 
but in the case that the exchange integral J has been made inhomogeneous^ that is, position dependent. The 
experimental means of achieving this J-inhomogeneity is left unspecified. One (perhaps naive) suggestion is 
that a sample of one antiferromagnetic material could be enriched, in places, with atoms from a different 
antiferromagnetic species. Note that since Jij > for all i, j, H is minimized when each spin anti-aligns with 
its nearest neighbours, just as in the homogeneous case. 

It has been known for many years that the dynamics of the homogeneous spin lattice (Jy = J > 0, 
constant) is described in the continuum limit by the 0(3) nonlinear sigma model. In it was shown, by 
adapting the dimerization process of Komincas and Papanicolaou 3^ , that this result generalizes in a natural 
way to the inhomogeneous problem considered here. That is, the continuum limit of system |^ is described 
by the nonlinear PDE 



where n : M^+i ^ ^3 y^^^ |^| ^ ^ 

is a rescaled time variable, and J{x, y) is a continuum approximant of Jij. 
The relationship between n and Sy- is subtle. If we choose a lattice spacing e > 0, assumed small, then the 
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spin dynamics associated with a solution n(t, x, y) of is 

Sy(T) = {-ly-^ sn{2V2seT,ie, je) + 0(e). (3) 

Roughly speaking, one should imagine the field n sampled on a square lattice of spacing e, partitioned into 
black and white sublattices, chess-board fashion. The white spins evolve as sn{t), the black as —sn{t), where 
t = 2V2seT. 

Equation Q is still the field equation of the nonlinear 0(3) sigma model on ]R^+^, but now spacctime is 
equipped with a spatially inhomogeneous Lorentzian metric, 

V = dt'-—^idx'+dy'). (4) 

Note that the static field equation is independent of J, so the usual Belavin-Polyakov lumps of the homogeneous 
system which we shall call bubbles in this context (in analogy with magnetic bubbles in ferromagnets), carry 
over unchanged to the inhomogeneous system. Given the relativistic nature of (jSJ, one can instantaneously 
boost these static bubbles to obtain moving bubble solutions. However, due to the inhomogeneity of rj, such 
bubbles do not necessarily travel along straight lines. A detailed study of the trajectory of a single bubble 
interacting with various J-inhomogeneities was presented in 0. The main conclusion is that, under most 
circumstances, bubble trajectories experience refraction, in analogy with geometric optics, J{x,y)~^ playing 
the role of the refractive index of the medium. The incident and exit angles of a bubble crossing a domain 
wall J-inhomogeneity, for example, are related by Snell's law, and total internal reflexion occurs if the impact 
is sufficiently oblique. 

In the current paper, we shall briefly review this key refraction phenomenon, then go on to study some 
new dynamical problems which were not considered in 0. Bubble refraction can be derived using only 
conservation laws and the assumption that the scattering of a bubble from a domain wall is approximately 
elastic 0- However, in this paper our discussion of bubble dynamics will be entirely within the gedoesic 
(or adiabatic) approximation of Manton 0, the standard approach to low-energy soliton dynamics in field 
theories of this type. The geodesic approximation, as applied to slow bubble dynamics, is set up in section |21 
This framework is then used in section |21 to study the interaction of a single bubble with domain wall, trough 
and disk J-inhomogeneities. Finally the rotationally equivariant dynamics of n coincident bubbles in a disk 
J-inhomogeneity is studied in sectional Here a dynamical phenomenon entirely absent from the homogeneous 
model, and with no optical analogue, is found. A charge n bubble set spinning about its symmetry axis can, in 
some circumstances, undergo periodic shape oscillations. These (approximate) solutions are spatially localized 
and, since the period of the shape oscillation need not be commensurate with the internal rotation period, in 
general, quasiperiodic in time. Hence they are quasibreathers. By contrast, a charge n bubble set spinning 
in the homogeneous model simply spreads out indefinitely. The inhomogeneity of J may thus act to stabilize 
spinning bubbles. 



2 The geodesic approximation 

Equation Q is the variational equation for the action 

5 = i y" dtdxdyy^\d^n ■ d^nrj^"" = J dt{T-V) 
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are identifed as the kinetic and potential energy functionals respectively. Note that while T depends on J, 
V does not. This is a consequence of the deep fact that the sigma model potential energy (often called the 
harmonic map or Dirichlet energy) on a space of dimension 2 depends only on the conformal class of the spatial 
metric. In our problem this metric is g = J{x, y)~'^{dx'^ + dy^) which is conformal to the Euclidean metric for 
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any choice of J{x,y). It follows that Belavin's and Polyakov's analysis of the static homogeneous model ^ 
carries over unchanged to the inhomogeneous system. There is a topological lower energy bound 

V[n]>47r\n\, (7) 

n being the topological degree of the map n : U {00} — > S*^. It is useful to define complex coordinates 
z — X + iy on the spatial plane and u = (rii + m2)/(l — n^) on the target sphere (the latter being the image 
of n under stereographic projection from (0,0, 1) to the equatorial plane). Then the bound for n > is 
attained if and only if u{z) is a rational map of algebraic degree n, that is 

ao + aiz-\ h a„z" 

u(z) = ; ; 8) 

where at, bi are complex constants, at least one of a„, 5„ is nonzero, and the numerator and denominator share 
no common roots. Since such maps minimize V globally within their homotopy class, they are automatically 
stable static solutions of the model. 

Without loss of generality we may choose the boundary value of n at spatial infinity to be (0, 0, 1), that is 
u{oo) — 00, so that 6„ = and hence a„ 7^ 0. Since a„ 7^ 0, we may divide the numerator and denominator 
by a„ and relabel a.i/a„ a^, bi/bn ^ bi, so that the general static solution is 

ao + aizH h a„_iz""i + z" 

^^(^) ~ — ; ; ; • (9) 

This is uniquely specified by the 2n complex constants ao, ■ ■ ■ ,fon-i- Hence the moduli space (or parameter 
space) of degree n static solutions, M„, may be identifed with an open subset of C^", namely the complement 
of the complex codimension 1 algebraic variety on which the numerator and denominator have common roots. 
For example, the general degree 1 bubble is 

uiz) = x~^e~"^{z-w), (10) 

X S (OjOo); ip G [0i27r] and w ~ wi + iw2 G C being constants interpreted as the bubble's width, internal 
phase and position respectively. Hence Mi is diffeomorphic to x C, where = C\{0}. In general M„ 
has real dimension 4n, and may be thought of as the configuration space of n unit bubbles. However, it is 
not diffeomorphic to (Mi)", and its parameters cannot globally be identified with the positions, widths and 
orientations of n separate bubbles. 

Following Ward (6| and Leese |7j , we may study the low energy dynamics of n bubbles within the geodesic 
approximation of Manton The idea is that M„, the 4n dimensional space of static degree n solutions, is 
the fiat valley bottom in the space of all degree n maps K.^ ^ S"^, on which V attains its minimum value of 
47rn. Consider the Cauchy problem in which the constituent bubbles of a static n-bubble are set off in relative 
motion with small speeds. Then the initial field lies in M„, the initial velocity field is tangential to M„, and 
the system has only a small amount of kinetic energy. Hence its subsequent motion in field configuration space 
is confined close to M„ by conservation of energy E = T + V (note that T is strictly positive). This suggests 
that a collective coordinate approximation wherein the motion is constrained to M„ for all time is sensible. 
This is often called an adiabatic approximation, since we assume that at each time the field configuration is 
well approximated by a static solution, but that the parameters in this solution may vary slowly with time. 

Let , . . . , 5*" be an arbitrary real coordinate system on M„ (for example, the real and imaginary parts of 
ao, . . . , a„_i, 60, ... , bn-i), and let us denote by uo{z; q) the static n-bubble corresponding to q = (q^, . . . , 9*"). 
Then the adiabatic approximation asserts that 

u{z,t) = uo{z;q{t)). (11) 

We now substitute (|ll|l into S, and obtain variational equations for q. Note that V = Airn, constant, for all 
fields of the form l(TT|l . and T is quadratic in time derivatives, so the action reduces to (Einstein summation 
convention applied) 

S^ljdt i^^^q^q^ - 4.n), 7.(.) - / il + \uoiz;qWyWW- ^''^ 
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This is the action for geodesic motion on M„ with respect to the metric 7 = "fij{q)dq^dq^ . Hence we expect 
u{z, t) to be well approximated by geodesic motion in (M„, 7), at least for low speeds. It should be emphasized 
that this does not mean that individual bubbles follow geodesies in physical space (M^,^), even in the case 
n = 1. Indeed, the motion of 2-bubbles even in Euclidean space (i.e. even in the homogeneous model) is highly 
nontrivial [HJITj- For example, two identical bubbles fired directly at one another do not pass through one 
another unchanged, or achieve some minimum separation then recede along their line of approach, as one might 
expect. Rather, they coalesce to form a ringlike coincident 2-bubble, then break apart and recede along a line 
perpendicular to their line of approach. This behaviour is generic to planar topological solitons in relativistic 
field theories, and is well-understood within the geodesic approximation 8 . The conceptual framework and 
validity of the geodesic approximation are discussed at length in 9 . 

The metric 7, and hence the dynamics, depends strongly on J{x,y). The main task in this approach to 
topological soliton dynamics is to calculate this metric. Even in the degree 1 case this is impossible to compute 
explicitly unless J takes a particularly simple, symmetric form. Nonetheless, much qualitative information 
about 7 can be deduced, and this can lead to a good understanding of the dynamics of a single bubble even 
in the absence of explicit formulae. In the case of higher degree, n > 2, the moduli space is too large for the 
general geodesic problem to be tractable, even if 7 is known explicitly. To make progress one must reduce the 
dimension by identifying totally geodesic submanifolds. This amounts to imposing symmetry constraints on 
the initial data. 

The metric in the case n = 1 has a particularly simple form. If J is constant, we are studying the standard 
sigma model, and it is found that x and ■0 are frozen by infinite inertia {■Jxx ~ ~ and 7 = AT:J~^dw dw, 
so bubbles just travel at constant velocity 0. At the other extreme, if J{x, y) — \ + + y"^, we effectively 
have the sigma model on a round two-sphere, and the bubble dynamics is much richer This choice has J 
unbounded, which is presumably unphysical, however. Assuming that J remains bounded, x and if) are frozen 
by an essentially identical argument to Ward's. Without loss of generality, we can set "ip — 0. The width x 
remains a free, but frozen, parameter. The metric on the width x, phase leaf of Mi, is then 

^ = Jx{w){dwl+dwl), (13) 

where 

. , , f dzdz 4x"^ f idudu 1 ^, 

/xH = / -TTT^ , .,-21, ..A,A2 = / , l.i2N2 t^.,.. , ...^2 ^ (14) 



J(z,z)2 (l-f X^^l-Z - Wp)2 J (1 + |-u|2)2 J(^u -f w)2 

the integral of J(x'u + w)^^ over the standard unit sphere on which u is a stereographic coordinate. Note that 

47r 

l™n^ " TT^^^"'! + '^"'2) 47rg. (15) 

In the limit of vanishing bubble width, therefore, the trajectories of isolated bubbles tend to geodesies, not 
just in M^, but also in the physical plane. One may regard the metric 7 on as a smeared out version of 
due to the bubble's finite core size. 



3 Bubble refraction 

The simplest J-inhomogeneity of all is the domain wall. Let us assume that J depends only on x, is constant 
outside a small neighbourhood (—6,6) of a: = and rises monotonically from J_ to J+ as x traverses this 
neighbourhood. Such a J{x) might arise from enriching one end of an antiferromagnet but not the other, for 
example. It follows immediately from H14() that the conformal factor fx{w) of the metric 7 on depends 
only on wi , is monotonically decreasing and has 

47r 

lini /x(wi) = -2". (16) 

tDl— »±00 

Our analysis of the interaction of the bubble with the domain wall will use only these properties of An 
explicit formula for can be obtained in the limit of vanishing domain wall width (5 — > 0. If we idealize the 
domain wall by a step function. 
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we find from l|14(l . and some elementary spherical geometry [21, that 



fli 



27r 



ai 



(18) 



Note that decreases smoothly from An/ J"^ to 47r/J^, as expected. 

To analyze geodesic fiow in M^, it is helpful to use the Hamiltonian formalism. The position coordinates 
wi,W2 have canonically conjugate momenta 



Pi = /(wi)wi, P2 = /(wi)w2, 



and the flow is generated by the Hamiltonian 



H = 



Viwi) 



(19) 



(20) 



Both H and p2 are integrals of the motion. Since geodesic trajectories are independent of initial speed, 
we may assume, without loss of generality, that H — ^. Trajectories now fall into a one-parameter family 
labelled by the conserved momentum p2, which given must satisfy p| < sup^^gR /(wi) = Att/J^. The 
image of a momentum p2 trajectory under projection to the (wi,pi) plane in phase space lies in the level set 
f{wi) — p\ = P2- A sketch of these (projected) trajectories is presented in figure^ They are confined to the 
bottle shaped region bounded by the curves pi = ±y/ f{wi) whose union is the level set with p"^ = 0. These 
curves themselves correspond to a bubble hitting the domain wall with incidence angle (i.e. orthogonal to 
the wall) and travelling straight through it, either left to right (upper curve) or right to left (lower). Two types 
of trajectory are evident. If < p| < Att/J^, then wi never vanishes and the bubble passes through the wall. 
We call these refracted trajectories. If An/J^ < p\ < A-k/J'^, the trajectories have \\mt^±ao wi{t) — — oo, so 
the bubble approaches the wall from the left and is reflected by it. We call these reflected trajectories. Let 
9{t) be the angle between the bubble's velocity and the x axis, so w — |u;|e'^(*^ Then 



sin0(i) = 



W2 



P2 



P2 



Hence for a refracted trajectory incident on the domain wall from the left, 



smf oo 



smf — oo 
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j: 



(21) 



(22) 



-oo) and 0{oo), are related by Snell's law of refraction, with 
oo) > J^l then p^ > An/J'\_ and total internal reflexion 



so the incident and exit angles of the trajectory, 6{ 
identifled with the refractive index. If sm9{— 
occurs. 

It is interesting to make a similar qualitative analysis in the case of a straight trough J inhomogeneity. 
Assume that J is again independent of y, that J{x) = Ji for \x\>l + 5, 5 > small; that J{x) = J2 < Ji for 
\x\ < 1 — <5; and that J interpolates monotonically between the constant values Ji, J2 in the narrow transition 
intervals (--1 — S, —1 + 6) and {1 — 6,1 + S). Given Snell's law for domain walls H22|l . one expects this to be a 
bubble guide: we have an infinite strip in which the refractive index is I/J2, greater than the refractive index 
1/ Ji in the surrounding region, so a bubble initially in and moving roughly along the strip should be trapped 
in the strip by total internal refiexion, just as a light ray is confined within a fibre-optic cable. 

The key point is, once more, to understand the conformal factor /^(w) of the metric on M^. Clearly 
depends only on Wi and tends to An/Jf as wi — > ±00. In general, it is hump shaped, with a single critical 
point, a global maximum. If we assume J(wi) is even, then so is f^{wi), and this maximum, which is always 
less than 47r/j|, occurs at wi — 0. Again, in the sharp wall limit, (5 — *■ 0, we can find an explicit expression 
for f^{wi). The point is that J{wi) is piecewise constant in this limit, so from ((T^ . we need only compute 
the area of the region on whose stereographic image is the vertical strip {u G C : — 1 < Re (xu + w) < 1}. 
This is easily deduced from the sharp domain wall calculation in (which led to formula (|18|l ^. yielding 



47r /I 1 
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(23) 
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Figure 1: Projected phase portrait of the geodesic flow in IVl!^ for a domain wall J-inhomogeneity. The exchange 
integral on the right of the wall, J+, exceeds that on the left, J_. Trajectories are bounded between the upper 
and lower bold curves, for which p2 ~ 0. Dashed trajectories, which have ^ < P2 < j?-, correspond to bubbles 
undergoing total internal reflexion. 



Note that this function has the quahtative features predicted. 

The Hamiltonian discussion of geodesic flow in (M^,7) is similar to the domain wall case. Once again we 
choose H = ^ and sketch the level curves fx{wi) — pi — P2 in the {wi,pi) plane for < p^ < /max = /x(0)- 
Again, two regimes emerge separated by the case P2 = pI — 47r/ J^, see figure|2 For p2 < p* the trajectories 
are unbounded and wi never vanishes. These correspond to a bubble hitting the trough from outside and 
passing right through it, being refracted twice in the process, so that its initial and final velocities coincide. 
For p2 > p* the (projected) trajectories are bounded and periodic. These correspond to a bubble being directed 
along the direction of the trough by total internal reflexion. The actual bubble trajectory in the physical plane 
is not periodic, of course, since it drifts uniformly in the y-direction {102 never vanishes; it is not constant, 
however). Note that even in this trapped regime, the bubble trajectory is not necessarily confined within the 
trough itself. Indeed, given any i? > 0, there exists a trapped geodesic for which maxj |wi(i)| > R. In this 
sense, the behaviour of bubble trajectories is quite different from that of light rays in geometric optics. 



P. 




Figure 2: Projected phase portrait of the geodesic flow in M^^ for a trough J-inhomogeneity: the exchange integral 
is Jl outside a vertical strip, within which it falls to J2 < Ji. The closed dashed trajectories, which have > 4?. 
correspond to bubbles being guided along the trip by total internal reflexion. 

The bubble guide effect does not depend on the trough being straight. The case of a circular trough (where 
the value of J is suppressed in an annulus) was considered in j2j, the results being qualitatively similar. 

Another interesting effect reminiscent of geometric optics is the phenomenon of bubble trajectory focussing 
by a circular "lens". Imagine J depends only on |z|, so that J(|z|) = J- for |z| < 1 — S, J{\z\) = J+ > for 
\z\ > 1 + 6, interpolating smoothly and monotonically between these values in the narrow annulus l~S < \z\ < 
1 + 6. We have a disk where the refractive index 1/ J_ exceeds that of the surrounding medium, so we expect 
single bubble trajectories to be refracted by this inhomogeneity like light rays passing through a circular lens. 
It follows immediately from H14(l that the conformal factor in the metric depends only on |w|. In the sharp 
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wall limit (5 ^ 0, it can be explicitly computed, 



1 - 



v/x* + 2x2(|u;|2 + l) + (|u;|2-l)2 



27r 



x^ + lwp-i 



(24) 

The geodesic problem in may easily be solved numerically, see figure 13 Note that parallel incident 
trajectories are approximately focussed by the lens, as expected. 




Figure 3: Focussing of parallel bubble trajectories by a disk J-inhomogeneity. The exchange integral is lower 
(J_ = 2) within the disk than in the surrounding plane (J_|_ = 3). 



4 Quasibreathers 

All the J-inhomogeneities we have considered enjoy some kind of symmetry, either translational (domain 
walls and troughs) or rotational (disks). In this section, we will continue to impose rotational invariance, so 
J = J(|z|). However, our emphasis will be very different: we will study rotationally equivariant n-bubble 
dynamics for n > 2. This corresponds to n coincident magnetic bubbles sitting, for all time, at the symmetry 
origin z = 0. This is dynamically interesting because for n > 2 the n-bubble width is no longer a frozen 
parameter. The n-bubble can spin internally and spread out or sharpen. So the dynamics of interest is no 
longer bubble transport, but rather the internal motion of the n-bubble. We will again model this within the 
geodesic approximation. 

First note that, if J = J(|z|), there is a natural isometric action of U{1) on (M„,7), given by 

u{z) >-> e-""M(e'"z). (25) 

The fixed point set of this action is necessarily a totally geodesic submanifold of (M„, 7), which we shall call 
{M'^,j\), where 7I is the induced metric. In fact fTJ, M'^ is, for all n, diffcomorphic to C^, since it consists 
of the rational maps 

m(z) = cz", ceC\ (26) 

Since U{1) equivariance is a valid symmetry reduction of the full field equation one may regard geodesic 
flow in {M'^j 7I) either as a low energy approximation of the symmetry-reduced field equation, or as a symmetry 
reduction of the general geodesic approximation. ^ 

The geometry of M^'^ is most easily understood by studying the lift 7 of 7I to the n-fold cover M"'^ of 
This is again diffcomorphic to C^, and the covering projection is 

n:Ml^^Ml\ n{p) = . (27) 

Hence, the n-bubble corresponding to p G is u{z) = {z/ p^^. Clearly any of the n points pe27rjfc/ri |y|cq^ 
/c = 0, 1, . . . , n — 1, give the same n-bubble. We may think of \p\ as the n-bubble width, and the argument 
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of p, modulo 27r/n, as being an internal phase. The potential energy density of an n-bubble is rotationally 
symmetric, vanishes at z = and \z\ — > oo, and attains its maximum at 

WHpI(^)''^ (28) 

The hfted metric on M^i is 

7 = /(H)dpdp, m) = 87rn'J^ (1 + r^^)^ ' ^^^^ 

In the homogeneous case, J{\z\) = Jq, constant, the metric is 7 = Cndpdp where C„ is a constant depending 
on n, so that {M'^,j) is actually isometric to the Euclidean punctured plane, and the lifted geodesies are 
simply straight lines traversed at constant speed. It follows that all geodesies have \p\ 00 either as i ^ 00 
or as t — > —00, and, except in the case where p{t) is a half ray hitting the missing point p = 0, both. The 
half rays correspond (when oriented inwards) to an n-bubble with no angular momentum collapsing to form a 
singularity in finite time. The straight lines missing p — correspond to n-bubbles with angular momentum 
which spread out indefinitely as |t| 00. In total, each spin n(z) in such a bubble executes n/2 complete 
rotations about the 713 axis. 

The behaviour is rather more interesting if we allow J(|2;|) to be non-constant. The "spherical" case 
J{\z\) = 1 + |zp was studied in detail in fl]. One nice, and rather unexpected, fact in this case is that the 
volume of M^^^, which turns out to be finite, is actually independent of n. One can isometrically embed 
in M.^ as a surface of revolution, and this surface has conical singularities of deficit angle 27r(l — n"^) at c = 
and c = 00. However, this case has J unbounded, which is not physically reasonable. We shall here consider 
the circular lens case, where J is suppressed in a disk centred on the origin. For simplicity, we will use the 
sharp wall limit, that is 

with J+ > We saw in the last section that the metric on could be computed explicitly for this J(|z|). 
We emphasize that this is not the metric of interest in the present context, however. Recall that position in 

is specified by the bubble's position in the physical plane, and the bubble width % is a frozen parameter. 
By contrast, position in M'^ is specified by the n-bubble's width and internal phase, and its location in the 
physical plane is fixed. In [2] an interesting bifurcation in the geodesic flow on was observed when x is 
varied: as x increases through a critical value, a pair of periodic geodesies appears in which the lump orbits 
in a circle centred on z = 0, one just inside the lens, one just outside. We will see that there is a geometrically 
similar bifurcation in the geodesic flow in M°^, but its physical meaning is now completely different, and the 
bifurcation parameter is not the bubble width (now a dynamical variable). 

Assuming J is given by H3U|I . one sees that 

27rn^ i TT ( ,T 



Since the integrand is rational, one can, in principle, write down / explicitly for any choice of n. The 
expressions are not terribly enlightening however, even for small n, so we will omit them. More useful is the 
fact that, provided J+/J- is not too large, M°i can be isometrically embedded as a surface of revolution in 
W" . Generating curves for these surfaces in the case n = 2, generated by the method described in JI], are 
depicted in figure 0] If J+/J_ = 1 the surface is just a (punctured) plane, swept out by rotating a straight 
half-line. As J+/J_ increases, the generating curve becomes hump-shaped, until eventually, it forms a bulge 
above a narrower neck. At this point, a bifurcation in the geodesic flow occurs. A pair of periodic geodesies 
with \p\ constant appears, and there are quasiperiodic (and periodic) geodesies which remain confined within 
the bulge. These correspond to quasibreathers. As J+/J- increases further, the isometric embedding is lost, 
but the existence of quasibreather geodesies persists. 



8 



-10 



-5 



° X ' 



10 



Figure 4: Generating curves for Mj'^ isometrically embedded as a surface of revolution in M.^ in the case of a sharp 
disk J-inhomogeneity with J^/J^ ~ 2, 3, . . . , 9 (top to bottom). The surfaces are swept out as these curves are 
rotated about the vertical axis a; = 0. Note that a neck develops as J+/J- becomes large. 

To confirm this visual reasoning, let us again consider the phase portrait of the geodesic flow. Introducing 
polar coordinates (cr, -0) on M'^ by ae^^ := p, the flow is generated by Hamiltonian 

the conjugate momenta being 

Pa = /(cr)o-, Pi, = (j'^f{a)ip. (33) 

Both H and are integrals of the motion, and we may choose H = ^ without loss of generality. The 
trajectories, after projection to the [(T,pa) half-plane (note cr > 0) are given by the level sets 

a\f{a)-pl)^pl (34) 

for G [0,oo). These level sets are plotted for J+/J- =2 and J+/J- = 4 in the case n = 2 in figure[Sl Note 
that as J+/J- exceeds a critical value (around 2.89 for n = 2), an island of periodic orbits appears in the 
phase portrait. This is due to the formation of a local maximum in (t^/((t) and hence, given the asymptotic 
behaviour ct^/(ct) ^ 47r J^^cr^ as tr — > oo, due to the formation of a pair of critical points of (T^/(cr). These 
points have a nice geometric interpretation, which we will now describe. First note that 

Now, the necessary and sufficient condition for M'^ to have an isometric immersion as a surface of revolution 
in is that —1 < S((t) < 1 for all a, see In this case, if we imagine the surface as being generated by 
rotating a generating curve in the xiX2 plane about the X3 axis (as depicted in figureEJ, then S(cr) = cos^(cr), 
where ^ is the angle between the tangent to the generating curve and the xi axis. Hence critical points of 
ct"^ f{a) are associated with points where the tangent to the generating curve is vertical, and a pair of such 
points develops precisely when the surface develops a pinched neck. This confirms our visual intuition that 
the island of closed trajectories appears when the pinched neck forms in M"'^. 
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Recall that the trajectories depicted in figure El are the projections of geodesies to the {a,pcr) plane. A 
closed trajectory in this picture does not necessarily correspond to a periodic geodesic: both the width a{t) and 
the phase e*'^*-*^ are periodic, but with (in general) incommensurate periods. These geodesies are therefore only 
quasiperiodic, so we call the oscillating n-bubbles associated with them quasibreathers. The time dependence 
of the width a and internal phase ?/; of a typical quasibreather is depicted in figure El Note that the shape 
of these solutions really does "breathe" . This is in contrast to the internally spinning magnetic solitons 
(sometimes called breathers) which occur in ferromagnetic spin lattices (with on-site anisotropy), where n^, 
has a fixed hump-shaped profile and the spins merely rotate at constant speed about the 713 axis jl2| . For 
n = 2, quasibreathers exist only for J+/J_ > 2.89 (approximately), but for any value of J+/J- > 1, there is 
some riniin > 2 such that the geodesic flow supports quasibreathers of any degree n > rimin- The dependence 
of rimin on J^/J^ is shown in figured 




Figure 6: The time dependence of the width a (solid curve) and internal phase (dashed curve) of a typical degree 
2 quasibreather, in a disk inhomogeneity with J+/J_ = 4. 
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Figure 7: The minimum allowed degree nmin of a quasibreather as a function of the inhomogeneity parameter 

J+/J_. 

5 Concluding remarks 

We have seen that, in the near continuum regime, a two-dimensional isotropic Heisenberg antiferromagnet 
with position-dependent exchange integral is described by the relativistic 0(3) sigma model on a spatially 
inhomogeneous spacetime. By conformal invariance, such a model supports static Belavin-Polyakov lumps, 
which we call bubbles. The trajectories of these bubbles are predicted to behave like light rays propagating in 
an inhomogeneous medium, with J^^, the inverse exchange integral, identified with the refractive index. In 
particular, bubbles incident on a straight domain wall are refracted in accordance with Snell's law, and total 
internal reflexion can occur. This refraction phenomenon was noted in 0. 

In this paper we have reviewed the analysis of bubble dynamics presented in [2], and extended it beyond the 
scattering problems previously considered to include the rotationally equivariant dynamics of higher degree 
bubbles. Here we found that disk J-inhomogeneities with sufRciently large J+/J_ support quasibreathers: 
degree n > 2 bubbles centred in the disk which spin internally while their shape oscillates with a generically 
incommensurate period. 

Our analysis has been conducted entirely within the framework of the geodesic approximation of Manton. 
The validity of this approximation for the 0(3) sigma model remains an open question. In the general context of 
relativistic field theories of Bogomol'nyi type, it has proved to be very succesful at describing soliton scattering 
processes at low to intermediate energies, where the solitons move essentially as free particles except during a 
relatively brief interaction phase 9 . The basic domain wall refraction effect seems likely to be a robust feature 
of the degree 1 bubble dynamics, therefore. 

The status of the quasibreather solutions is less straightforward. It is unlikely that the sigma model will 
support genuine quasiperiodic solutions for any bounded choice of J(|z|), because one expects an oscillating, 
spinning bubble to radiate energy to infinity in the form of small amplitude travelling waves. So the quasiperi- 
odic geodesies in M'^^ found here cannot be globally close to the true solutions of the model with the same 
initial data. One should think of the geodesic approximation as predicting that long-lived oscillating, spinning 
bubbles exist in the model which are approximately quasiperiodic, but decay over a long time scale. 

The only geodesies for which rigorous error estimates have been found are the radial half rays p(t) = —vt 
in M^'^ for n > 3 (in the homogeneous case, J constant, though this is unlikely to be crucial). These geodesies 
describe an equivariant n-bubble collapsing in finite time. Rodnianski and Sterbenz |13| have shown that 
the actual solutions do indeed collapse in finite time, but that the collapse rate predicted by the geodesic 
approximation receives a logarithmic correction which becomes significant when the bubble is extremely close 
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to collapse (i.e. extremely narrow). In the present context, the significance of this deviation is debatable: when 
the bubble is very narrow, the continuum approximation which led to the sigma model in the first place has 
presumably already broken down. Indeed, the most important test of the dynamical phenomena predicted 
here and in 1] is not provided by rigorous analysis of the PDE ^ , but by direct numerical simulation of the 
spin lattice itself. Such simulation is computationally far more intensive than the methods we have used 
here, and lies beyond the scope of this paper. 
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